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Search for Glueballs from Three-body Annihilation of pp in-Flight 

Bing-Song Zou a * 

a Queen Mary and Westfield College, London El 4NS, UK 

Lattice QCD and other theoretical models predict that the 0~ + , 2 ++ and 2 h glueballs 
have masses in the range of 2.0 to 2.4 GeV. For resonances in such an energy range, three- 
body decay modes are expected to be large. The strategy for looking for these glueballs 
and the newest results from studying Crystal Barrel data on the three-body annihilation 
of pp in flight are presented. 

1. STRATEGY FOR GLUEBALL HUNTING IN pp ANNIHILATION 

Proton-antiproton annihilation is regarded as a favorable process for glueball produc- 
tion. For pp annihilation, there are two possible ways to produce glueballs as shown in 
Fig.l(a&b), which are so called "production" and "formation" mechanisms, respectively. 



P 





(o) 
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Figure 1. Glueball production mechanisms for pp annihilation: (a) "production" mech- 
anism; (b) "formation" mechanism. Dashed and helix lines are for mesons and gluons, 
respectively. 



For the ++ glueball ground state, lattice QCD predicts its mass to be 1.45 ~ 1.8GeV/c 2 
||. For a glueball in such a mass range to be produced from pp annihilation, it can only 
come from the "production" mechanism. Indeed, by studying pp — > 37r° & ir°r]r], Crystal 
Barrel Collaboration discovered the /o(1500) resonance which is now regarded as the 
best ++ glueball candidate 0. There is also some new evidence for / (1770)|{||]. 

If the /o(1500) is really a glueball, it sets a mass scale for glueballs of other quantum 
numbers. The 0~ + , 2 ++ and 2~ + glueballs are predicted to be around 2.1 ~ 2.4 GeV by 
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various theoretical models |l],|5|||. For resonances in such an energy range, they should 
be mainly produced from the "formation" mechanism and three-body decay modes are 
expected to be large. 

Crystal Barrel has taken a lot of data for all neutral final states in flight. Then the 
question is which three body channels we should study first. Here we can take some 
lessons from charmonium decays. For rj c (0 h ), Xco(0 ++ ), Xc2(2 ++ ) hadronic decays and 
J/ty radiative decays, they definitively go through two-gluon intermediate states. The r), 
rf, a and / (1500) seem to be favoured decay products of two-gluon states |7|||. The r]a, 
rj'a and 77/0 (1500) are expected to be large decay modes of 0~ + glueballs while the 77/2 
and t//2 are expected to be large for 2~ + and 2 ++ glueball decays. These decay modes 
have 7r°7r°77, ir ir T]' and 3rj as their final states. 

Hence our strategy for hunting _+ , 2 ++ and 2~ + glueballs is to study resonances formed 
by pp and decaying into 7r°7r°77, 7T°TT°r]' and 3r] final states first. 

2. STATUS OF THREE-BODY ANNIHILATION IN-FLIGHT 

Crystal Barrel at LEAR has collected data triggering on neutral final states at beam 
momenta 0.6, 0.9, 1.05, 1.2, 1.35, 1.525, 1.642, 1.8 and 1.94 GeV/c, which correspond to 
center-of-mass energies ranging from 1.96 to 2.41 GeV/c 2 . An average of 8.5 million all 
neutral events were taken at each momentum. 

These data have been processed. Rough number of selected events at each beam mo- 
mentum and background level for reconstructed three-body channels from 67 and 77 
events are listed in Table 1. Among them, the background level for the tt tc uj channel 
has not been investigated yet; the background level for 7t n°ri' is too high to do partial 
wave analysis. So we have been analyzing the first four channels. 

Table 1 

Rough number of selected events at each beam momentum and background level for 
reconstructed three-body channels from 67 and 77 events. 



channel 


Number of Events 


background level 


7T U 7T U 7r U 


~ 150A 


1% 




~ 70A 


3% 


7T°r7?7 


~ QK 


6% 


i]i]ri 


~ 150 


3% 


Ay 


~ IK 


50% 




~ 70A 





For 7r 7r°7r° channel, the most obvious contributions come from /2(1270)r/ and fo(1500)i] 
intermediate states. The partial wave analysis is in progress. 

For 7r°rjrj channel, the statistics is not enough for a full partial analysis including both 
production and formation amplitudes. Some effective formalism was used to concen- 
trate on searching for resonances in the production mechanism[|J. The main results are: 
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/o(1500) -> 7777 is clearly seen; / (1750 ~ 1800), / (2100) and a 2 (1660) are confirmed; 
in addition there is a broad / 2 (1980) — > 7777 with mass M = 1980 ± 50 MeV and width 
r = 500 ± 100 MeV. 

For 777777 channel, the statistics is low. But we can still learn something. In Fig.2, we 
show the real data and some Monte Carlo Dalitz plots for pp — > 777777 at 1.8 GeV/c. 
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Figure 2. Real data and Monte Carlo Dalitz plots for pp — >■ 777777 at 1.8 GeV/c. 



From the real data Dalitz plot in Fig.2, the / (1500) is obviously there. For the beam 
momentum of 1.8 GeV/c, the orbital angular momentum between / (1500) and 77 is ex- 
pected to be < 2, which corresponds to the initial states of 0~(Z/ = 0), l + (// = 1) and 
2^ (If = 2). From their Monte Carlo Dalitz plots in Fig.2, without fitting data, it is 
already clear that the 0~ — > /o(1500)77 is the most obvious contribution to pp — > 777777. 

Preliminary results of the cross section for pp — > 777777 is shown in Fig. 3. There are two 
peaks at about 2.15 GeV and 2.33 GeV. They may be due to statistical fluctuations. But 
if they are due to two resonances, they are most likely 0~ + resonances with the / (1500)?7 
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decay mode. 




The only channel which is favourable for hunting 0~ + , 2 ++ and 2~ + glueballs and has 
enough statistics for a full amplitude partial wave analysis is the 7r 7r°?7 channel. Its 
cross section is shown in Fig.4. There are clear enhancements at around 2.05 and 2.3 
GeV. Note that for a constant amplitude the cross section should decrease steadily as the 
energy increases. 

For the channel, projections on to M(tcit) and M(7ttj) at 900, 1200, 1525 and 

1800 MeV/c are shown in Fig.5. The / 2 (1270), a (980) and a 2 (1320) are clearly visible. 
The /o(980) and /o(1500) are visible on the M(mi) projection, but rather weak. As beam 
momentum increases, the / 2 (1270) peak becomes stronger while the a 2 (1320) peak gets 
weaker; this is a natural reflection of the rapidly opening phase space for / 2 (1270)?y, whose 
threshold is at a mass of 1820 MeV. 

Based on these data, a full amplitude analysis describing both production and decay of 
these resonances is carried out for each moment um||. 



3. AMPLITUDE ANALYSIS OF pp -> n°n°r] 

For the 7t 7t i] final state, possible pp initial states are _+ , 2~ + , 4 _+ etc. for pp spin 
singlet states, and 1 ++ , 2 ++ , 3 ++ , 4 ++ , 5 ++ etc. for pp spin triplet states. For our case 
with center-of-mass energies below 2.41 GeV, only 0~ + , 2~ + , 1 ++ , 2 ++ , 3 ++ and 4 ++ are 



expected to be significant \W\ and this has been confirmed in our analysis; 4 has been 
tried, but is not significant. Their corresponding pp total angular momentum J, orbital 
angular momentum L and total spin angular momentum S in the usual contracted form 
2S+1 Lj are: % for 0~+ 1 D 2 for 2~+, 3 Py for 1++, 3 P 2 or 3 F 2 for 2++, 3 F 3 for 3++, and 
3 F 4 or 3 H 4 for 4++. 

Let us choose the reaction rest frame with the z axis along the p beam direction. Then 
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Figure 5. Data (points with error bars) and fit (solid line) of invariant mass spectra for 
7r°7T (1 entry/event) and 7i°r] (2 entries/event). Beam momenta are given in each panel. 
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the squared modulus of the total transition amplitude is the following ||1 If: 

I = \A - + + A 2 - + \ 2 + {A^ 1 + A^ + t\ 2 + {A^T 1 + A^=~ li2 



j_ I aM=0 i aM=Q\2 i I aM=1 i aM=\\2 , I aM=-1 , /|M=-1|2 /i\ 
~H /i 2++ ~+~ ^4++ | ~r 1^2++ ~+~ ^4++ | + | /i 2++ ~r ^4++ I \ l I 

+2Re[(A^ + AffiiAg? + Af + tY - (^T 1 + AjJr 1 )^" 1 + A^" 1 )*] 

where M is the spin projection on the z-axis. Partial wave amplitudes Ajpc are con- 
structed from relativistic Lorentz covariant tensors, Breit-Wigner functions and Blatt- 
Weisskopf barrier factors flj"2"| , |i"5| . The barrier factors use a radius of 1 fm. The ^(1270)77, 
a2(1320)7r, ao(980)7r, arj, fo(980)rj and /o(1500)^ intermediate states are considered. The 
/o(980) is fitted with a Flatte formula using parameters determined previously |l4j]. The a 



is fitted with the parameterization A of Ref. [|14]]. Other resonances are fitted with simple 
Breit-Wigner amplitudes using constant widths. Full formulae and additional details will 
be given in Ref.[Q. Based on these formulae, the data at each momentum are fitted by 
the maximum likelihood method. 

The fit is shown in Fig. 5 for the mass spectra for beam momenta at 900, 1200, 1525 
and 1800 MeV/c. The quality of the fit for other beam momenta is similar. The fit to 
projections is obviously not perfect. These may be due to additional components such 
as ao(1450)7r, a 2 (1660)7r and even p(1405)7r intermediate states. An angular momentum 
decomposition of this small effect is not possible, but fits including such states produce 
little effect on the dominant components. Since we are mainly interested in scanning the 
larger components from / 2 (1270)?y, arj, a 2 ~K and ao(980)7r intermediate states, we ignore 
those smaller contributions for the present study. 

The data points with error bars shown in Fig. 6 are our fitted results for the partial wave 
cross sections at each momentum for pp — > 7i -n-°i] with rj 77. Only those partial waves 
with significant contributions are presented. There are rich peak and dip structures. 

For 4 ++ , a peak around 2100 MeV is clear for all 4 ++ partial waves. It is fitted by 
a Breit-Wigner amplitude with the mass and width fixed to the PDG values for the 
well established 4 + resonance ^4(2050) [0. The shift of the peak position is due to the 
centrifugal barrier factors for both initial and final states. It appears in / 2 r7 and a 2 7r 
with comparable strength. In addition to the /4(2050), there is clearly another 4 ++ peak 
around 2320 MeV in 4 + — > f^r] in the M=l partial wave. This may be identified with the 
/ 4 (2300) listed in the Particle Data Tables 0. 

A 4 ++ resonance around this mass has also been observed be VES in ?77r + 7r _ in the nA 
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For 2 ++ , two peaks around 2020 MeV and 2350 MeV have masses and widths compatible 
with / 2 (2010) and / 2 (2340) listed by the PDG [|7| as established particles. In addition, a 
peak around 2230 MeV shows up clearly in the / 2 7/ mode. It has a mass compatible with 
£(2230) observed in J/\l/ radiative decays |7|Jl^l, but has a larger width of ~ 150 MeV. 



For 2 _+ and 3 ++ , both have two peaks around 2050 and 2300 MeV. No corresponding 
entries exist for them as yet in the Particle Data Tables 0. For 1 ++ , there is a strong 
enhancement at the low energy end, decaying dominantly into a 2 7r. 

For h , there seems to be a broad component plus a peak in rja at ~ 2140 MeV with 
width ~ 150 MeV. The broad component may correspond to the broad 0~ + object used 
in describing J/ty radiative decays to pp, luu, K*K*, (jxp and rj7m |T7|]. The peak at 2140 
MeV decays dominantly into arj. 
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Figure 6. Cross sections for partial waves of significant contributions to pp — > 7r°7r°r/ with 
7] — > 77. For diagrams with two components, the first label (full line) corresponds to the 
bigger component. The curves are the fit of Breit-Wigner amplitudes to the data points 
in the figure and the relative phases between components. 
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Table 2 

Summary of fitted masses, widths and branching ratios corrected for their unseen decay 
modes. The mass and width of /4(2050) are fixed at PDG values, and the status of the 
CT state at 2140 MeV is questionable, as discussed in the text. The /i(1700) is beyond 
the accessible mass range. All states have I = 0, G = +1. 



j PC 


M(MeV) 


r(Mev) 


^ pp^ fov i n3 

1 tot 


wi 2 10 J 

1 tnt 


% ^10 

1 tot 


10 

1 to* 


4++ 


2044 


208 


0.54 ±0.14 


5.1 ±0.8 


- 


- 


4++ 


2320 ± 30 


220 ± 30 


1.3 ±0.4 


1.0 ± 1.0 


- 


- 


_| L 

3 ++ 


2000 ± 40 


250 ± 40 


0.12 ± 0.08 


0.6 ± 0.6 




0.23 ± 0.11 


3++ 


2280 ± 30 


210 ±30 


1.7 ±0.4 


4.5 ±2.6 




0.23 ±0.19 


2++ 


2020 ± 50 


200 ± 70 


2.1 ±0.4 


4.3 ± 1.2 






2++ 


2240 ± 40 


170 ± 50 


2.5 ±0.6 


1.6 ± 1.6 






2++ 


2370 ± 50 


320 ± 50 


0.88 ±0.64 


16 ±5 






1++ 


~ 1700 


~ 270 










1++ 


2340 ± 40 


340 ± 40 


0.6 ±0.6 


60 ±30 




0.84 ±0.53 


o-+ 


2140 ±30 


150 ± 30 


1.9 ±1.7 


6.0 ±6.0 


10 ±5 




2-+ 


2040 ± 40 


190 ± 40 


3.0 ±0.3 


5.0 ±2.1 




0.4 ±0.2 


2-+ 


2300 ± 40 


270 ± 40 


2.8 ±0.7 


2.0 ±2.0 




0.5 ±0.5 



Using interfering sums of the Breit-Wigner amplitudes to fit the partial wave cross 
sections in Fig. 6 as well as relative phases between different partial waves, we obtain 
the masses, widths and branching ratios shown in Table 2. The corresponding Argand 
plots for all partial waves are shown in Fig. 7. Besides the obvious resonances mentioned 
in previous paragraphs, we need another 1 ++ resonance at about 2340 MeV with width 
~ 270 MeV. Without it, we cannot describe the relative phase between 1 ++ and 4 ++ 
partial waves; also we would need the lower 1 ++ resonance to be very narrow (< 50 MeV) 
in order to explain the sharp increase in the 1 ++ partial wave cross section at low mass. In 
our present fit with two 1 ++ resonances, the /i(2340) amplitude interferes destructively 
with the tail of the lower 1 ++ resonance and causes the sharply decreasing cross section 
with a broad dip around 2340 MeV. The phase motion caused by this /i(2340) can be 
seen clearly in the Argand plot for the 1 ++ — > a 2 n partial wave in Fig. 7. 

In Table 2, the branching ratios are calculated at the resonance masses and are corrected 
for their unseen decay modes, except for ao(980) where T aQ7T = r ao ^ vn7r . 

Among these resonances, the 7/(2140) and ^(2240) look special. Both have relative 
narrow decay width. The ?7(2140) decays dominantly into r/a; the / 2 (2240) has the largest 
f2V/ a 2ft ratio. These properties suggest that they may have larger mixing of glue com- 
ponents than other resonances. 

4. SUMMARY AND OUTLOOK 

In summary, from a full amplitude analysis of pp — > 7r°7r°?7 in-flight, we have observed 
a new decay mode riirir for three established resonances / 4 (2050), / 2 (2010) and / 2 (2340). 
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In addition, we have observed 8 new or poorly established resonances in the energy range 
from 1960 to 2410 MeV, i.e., / 4 (2320), / 3 (2000), / 3 (2280), / 2 (2240), A (2340), 77(2140), 
^(2040) and ^2(2300). Among them, the h 77(2150) has very large decay branching 
ratio to rja; 2 ++ / 2 (2230) has the largest f2i]/ a 2^ ratio; both have relative narrow total 
decay width. These properties suggest that they may have larger glueball components. 

For a further study of the 77(2140) and / 2 (2240), we are going to reconstruct ir°iT ri' 
channel from IO7 events, where we expect less contamination from other channels. The 
main purpose here is to scan the f 2 r]' and ar]' modes which are also expected to be 
favourable decay modes of gluballs. 

We are also going to scan / 2 7r° from 7r°7r 7r° final state to study isovector qq states and 
scan f 2 uj from 7r°7r°cj final state to study isoscalar qq states. Both 7r°7r 7r° and 7r°7r°c<j 
channels have enough statistics for a full amplitude partial wave analysis and have f 2 
band as their most obvious contribution in their Dalitz plots. 

From these analyses of three-body annihilation in-flight, combined with information 
from two-body annihilations, we hope to establish the 2.0 ~ 2AGeV/c 2 meson spec- 
troscopy which is crucial for identifying gluballs and understanding quark confinement. 

Acknowledgement: I thank the organizers of the LEAP98 conference for the invi- 
tation to give the talk in such a nice place, Sardinia. I am very grateful to D.V.Bugg 
and A.V.Sarantsev for the fruitful collaboration on data analysis, and to Crystal Barrel 
Collaboration for producing the beautiful data. 
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